Polariton condensation and lasing in optical microcavities - the decoherence driven 

crossover 

M. H. Szymanska, P. B. Littlewood, and B. D. Simons 
Theory of Condensed Matter, Cavendish Laboratory, Cambridge CB3 OHE, UK 

(February 6, 2008) 

We explore the behaviour of a system which consists of a photon mode dipole coupled to a 
medium of two-level oscillators in a microcavity in the presence of decoherence. We consider two 
types of decoherence processes which are analogous to magnetic and non-magnetic impurities in 
superconductors. We study different phases of this system as the decoherence strength and the exci- 
tation density is changed. For a low decoherence we obtain a polariton condensate with comparable 
excitonic and photonic parts at low densities and a BCS-like state with bigger photon component 
due to the fermionic phase space filling effect at high densities. In both cases there is a large gap in 
the density of states. As the decoherence is increased the gap is broadened and suppressed, resulting 
in a gapless condensate and finally a suppression of the coherence in a low density regime and a 
laser at high density limit. A crossover between these regimes is studied in a self-consistent way 
analogous to the Abrikosov and Gor'kov theory of gapless superconductivity [1]. 
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I. INTRODUCTION 

The miniaturisation and improvement in the quality 
of optical cavities in recent years led to the achievement 
of strong-coupling regime of light-matter interaction in 
many physical systems. The strong coupling regime is 
characterised by well-developed coupled modes of light 
and electronic excitations, called polaritons. Polariton 
splitting has been experimentally observed for atoms [2] , 
quantum wells [3] and bulk excitons [4], excitons in or- 
ganic semiconductors [5,6], exciton complexes [7] and 
glass spheres. The strong coupling regime has also been 
achieved for coupled Josephson junctions [8] in a mi- 
crowave cavity. 

With well- developed modes, that are sharp and have 
long lifetimes, a natural question becomes the existence 
of coherent, condensed states. A theoretically con- 
structed polariton condensate is a mixture of coherent 
state of light and coherent state of massive particles in 
the media. It is characterised by two order parameters: 
the coherent polarisation and the coherent photon field 
and exhibit a gap in the excitation spectrum [9]. Since 
the polariton condensate would be a source of coherent 
light the natural question arises how it is different from 
and how it can be connected to the traditional laser. 

The laser is a weak-coupling phenomenon: a coherent 
state of photons created by stimulated emission from an 
inverted electronic population due to strong pumping. 
The polarisation of the medium is heavily damped and 
the atomic coherence is practically zero. A coherent pho- 
ton field, oscillating at the bare cavity mode frequency, 
is the only order parameter in the system [10]. 

The crossover between a polariton condensate and a 
laser is sometimes mistakenly attributed to an increase 
in density and a crossover from bosonic (exciton) degrees 
of freedom to fermionic particle- hole pairs. The absence 



of polariton splitting is associated with disappearance of 
coherence, which is why experiments are concentrated at 
low densities. In this work we show that the trend is 
an opposite one - the condensate is more robust at high 
densities. One needs to remember that polaritons, and 
so polariton splitting, are normal state excitation and so 
disappearance of polariton splitting in the normal state 
does not indicate what would happen if the system was 
condensed. 

The issue of electronic coherence in polariton con- 
densate is independent of whether the excitations are 
more "excitonic" or more like a two-component plasma 
(ionised electron-hole pairs), or indeed whether or not 
there is saturation (phase-space filling) of the electronic 
states. It has been shown [9] that although the increase 
in the density of electronic excitation leads to nonlinear- 
ities in the polariton system which cause the collapse of 
the splitting between the two polariton peaks in the nor- 
mal state it does not destroy condensation even at very 
high excitation densities. The saturation in the fermionic 
space forces the condensate to become more photon like 
as the excitation density is increased and to have more 
BCS-like character, but nevertheless the coherence in the 
media and the gap in the excitation spectrum are present. 
The change in the density of the electronic excitations 
leads to the crossover between the regime of BEC of po- 
laritons and collective BCS-like state 

This article will argue that the real enemy of conden- 
sation is decoherence, not density, and that it is pre- 
cisely decoherence which drives the polariton conden- 
sate towards the laser regime. We show that only the 
self-consistent inclusion of decoherence processes allows 
to establish a crossover between an isolated condensate 
and a laser. The widely used quantum Maxwell-Bloch 
(Langevin) equations with constant decay rate for po- 
larisation are not correct in a regime when the coherent 
polarisation is large and the gap in the density of states is 
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present. We develop a self-consistent method analogous 
to the Abrikosov and Gor'kov theory of gapless super- 
conductors, which allows us to study the stability of the 
polariton condensate at low decoherence strength and to 
established the crossover to laser behaviour as the deco- 
herence is increased. 



II. MODEL 

The model we consider in this work is schematically 
shown in Figure 1. It consists of a set of N two-level oscil- 



lators dipole coupled to a single mode of electromagnetic 
field confined in a cavity. This system is then subject to 
various decoherence, pumping and damping processes. 
These processes can be of a different physical nature, de- 
pending on the material, but their exact details are not 
that important for a general theory. They can be de- 
scribed, similarly to laser theory, as baths of harmonic 
oscillators coupled to the system in a way that gives the 
same effect as the real physical environment. 
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FIG. 1. Sketch of the model studied in this work: the system of two-level oscillators dipole coupled to a single cavity mode 
interacting with various types of environment. 



Our model takes into account the major Coulomb in- 
teraction between the electron and hole within the ex- 
citon, the phase-space filling effect, disorder in the ma- 
terial (inhomogeneous broadening of excitonic energies) 
and various types of decoherence effects. However, it 
does not include screening and Coulomb interactions be- 
tween excitons. Therefore, it gives a very good descrip- 
tion of tightly bound, Frenkel - type of excitons localised 
by disorder or bound on impurities, molecular excitons 
in organic materials, atoms in the solid state or Joseph- 
son junctions arrays in microwave microcavity. And it 
gives only a qualitative description within a mean-field 
approximation for other types of excitons like Wannier 
excitons or excitons propagating in a sample (see Section 
VIII). 

The model we use is the minimal required to describe 



the essential physics. At this stage we do not intend 
to model any particular medium with its complex inter- 
actions which would only make the general picture less 
clear. More specific details of the particular medium 
could, however, be included straightforwardly into the 
formalism. We will discuss this possibility in the Sec- 
tion VIII. 



In the following, we consider a system comprised of 
an ensemble of N two-level oscillators with an energy ej 
dipole coupled to a single cavity mode. The correspond- 
ing microscopic Hamiltonian takes the form 
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H s = uj c ^tp + ^ ej(b'jbj - a]aj) 
3=1 

N 

+ 2-^(6^ + ^.) (1) 

3 = 1 

where Fermionic operators bj and a d - annihilate electrons 
in the upper and lower states respectively, while the 
Bosonic operator ip annihilates the photon. Here the sum 
extends over the possible sites j where an exciton can be 
present (e.g. different molecules or localised states associ- 
ated with a disorder potential). Matrix elements gj/VN 
describe the interaction of the photon with the two-level 
oscillators. 

Generally, the effect of the environment on the be- 
haviour of the system can be modelled through the in- 
teraction of the internal degrees of freedom with a bath. 
Taking into account different physical processes, the most 
general coupling is of the form 

Hsb = J29k(^d k + d{?P) 
k 

+ E [ b j a j(9jk<i,k + 9- k c-,k) + h.c 
jk 

+ J2 T ?k(fy j +a]a j )(cl k +c 1 , k ) 

jk 

+E r S ) ( 6 5 6 i-°5 o i)(4 lfc +c2,fc). ( 2 ) 

jk 

where 

k 

+Ul,kc\ tk Ci ik +W2,fc4 |fe C 2 ,fc 

describes the Hamiltonian of the bath. Here the differ- 
ent modes k of the bath are indexed by (independent) 
Bosonic field operators dk, c+,fc, c_ ; fc, ci,fc and C2,fe. Here 
the first term in (2) describes the decay of the photon field 
from the cavity. Matrix elements g^ k describe the inco- 
herent pumping of two-level oscillators, while the matrix 
elements g^ k contain all of the higher energy processes 
which destroy the electronic excitations such as the ra- 
diative decay into photon modes different from the cav- 
ity mode. Apart from their dephasing effect, together, 
these processes cause a flow of energy through the sys- 
tem. However, in the steady state, the total number of 
excitations 

3 

the sum of photons and excited two-level oscillators, is 
constant. Finally, the third and fourth terms describe all 
those lower energy dephasing processes, such as collisions 



and interactions with phonons and impurities, which con- 
serve the total number of excitations in the cavity. Such 
contributions can be divided into a part which act sym- 
metrically (r^) or antisymmetrically (r^ 2 )) on the upper 
and lower levels. Altogether these four terms contain all 
the essential mechanisms of decoherence. 

To assimilate the effect of the different mechanisms 
of decoherence, we will find it useful both intuitivelty 
as well as techincally (see later) to draw on an analogy 
between the Hamiltonian of the system and that of a su- 
perconductor [14,15]. Referring to the states indexed by 
the Fermionic operators bj as 'particle-like', and those in- 
dexed by Oj as 'hole- like', Hs can be interpreted as a BCS 
Hamiltonian for a superconductor with an imposed ho- 
mogeneous 'superconducting' order parameter ip. With 
this analogy, it is clear that the second and third terms 
of Hsb (2) affect a mechanism of 'pair-breaking', while 
the fourth term is compatibile with the symmetry of the 
Cooper pairs. The former act on the system as dynami- 
cally fluctuating magnetic impurities while the latter de- 
scribe dynamical fluctuations of a normal non-magnetic 
potential providing only an inhomogeneous broadening 
of the energies. 

When the pumping and photon decay rates are high, 
the system would be driven out of equilibrium. However, 
if the thermalisation rate is in excess of the speed at 
which the system is pumped, an equilibrium assumption 
can be justified. In this work we will limit our considera- 
tions to this regime focussing on the effect of decoherence 
on the equilibrium system. Choosing both the pumping 
and decay rates to be small — allowing thermal equilibra- 
tion their ratio can be used to fix the total excitation 
density n ex . Nevertheless, within this quasi-equilibrium 
regime, both the third and fourth terms in (2) can, in 
fact, be arbitrary large since they do not couple to the 
total density n ex . There is no restriction on the density 
of excitations nor on the decoherence rate. 

Therefore, on this background, we will consider the 
total Hamiltonian 

H = H S + H SB +H B , (4) 

where 

Hsb = Y, T %\b]b 3 + a) aj )(4 tk + c hk ) 

jk 

+ E r S } vfa a ) a M.k + <*,*). (5) 

jk 

includes only the includes only the decoherence mecha- 
nisms which conserve n ex , while 

r b = E E^< fcC ifc Clfe - 

i=l,2 k 

Although, in general, the coupling constants can be 
site dependent, for simplicity, we will suppose that the 
modes of the bath couple with equal strength to the two- 

(i) (i) 

level systems setting T\ k ' i — > T^. . Similarly, in the fol- 
lowing, we will assume that the coupling of the two-level 
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systems to the cavity photon is independent of the site 
index j, g 3 ^ g. 

When the interactions between the environment and 
the system are large, the Hamiltonian above provides 
the basis of the standard theory of lasers (see, e.g., 
Refs. [10,11]). At the same time, if we set the coupling 
constants between the system and the environment to 
zero, the ground state of the Hamiltonian H$ forms a 
polariton condensate. Thus by varying the magnitude 
of the coupling between the system and the baths the 
model provides the means to move smoothly between an 
isolated condensate and other phases driven by the deco- 
herence. This provides a means to explore the stability of 
the polariton condensate to interactions with the outside 
world at a small coupling strength, and to establish the 
connection between polariton condensation and lasers as 
the decoherence is increased. 

The standard assumption (which in many cases is 
physically correct) is that the environment leads to rapid 
dephasing of the exciton polarisation (i.e. Ti is very 
short). It is thus generally assumed that the polarisation 
is very small, and that the coherent photon field is the 
dominant order parameter. In such a case the baths can 
be averaged out before the exciton-photon interaction is 
studied. This in turn leads to the well-known quantum 
Maxwell-Bloch (Langevin) equations for the photon field 
and polarisation, essentially of the form 

± { a^)=i([Hs^b])-^b) . (6) 

Crucially, Eq. (6) can be derived by assuming that the 
lifetime T 2 for polarisation of a single two-level oscillator 
is the same as that of the macroscopic ensemble of two- 
level oscillators (which would correspond to introducing 
a separate bath for each two- level system). However, in 
general Eq. (6) is not correct [12]. When T 2 is long, the 
macroscopic ensemble of two-level oscillators can exist in 
a collective state characterised by a large coherent po- 
larisation, and the assumptions which lead to (6) cannot 
be justified. Moreover, the constant decay rate 1/T 2 in 
(6) is a critical parameter; even at arbitrarily small deco- 
herence it leads to completely different solutions from 
those in the absence of an environment (as shown in 
Rcf. [12]). This criticality is however unphysical and 
arises only due to approximations used in deriving the 
Langevin equations. In fact, this conclusion seems not 
to be widely appreciated. Indeed, in a relatively recent 
publication [13] similar decay constants in the equations 
of motion for the coherently driven excitonic insulator 
have been used, leading to the conclusion that the exci- 
tonic insulator phase cannot exist for an arbitrary small 
decoherence. 

One can gain some physical insight into this problem 
by considering the evolution of the density of states: The 
ideal condensate has a gap in the density of states which 
would still be present for small decoherence. It is evident 
that the coherent fields in this regime cannot be damped 



just by constant decay rates independent of frequency as 
there are no available states in which to decay. As the 
decoherence is increased this gap gets smaller and finally 
is completely suppressed causing the coherent fields to be 
strongly damped, as in lasers. In this regime Eq. (6) is 
perfectly valid. However to be able to study a crossover 
from the fully phase coherent polariton condensate to 
a laser one needs to include the environment in a self- 
consistent way in which the possible gap in excitation 
spectrum and a large polarisation would be taken into 
account. 

The analogy between superconductivity and the exci- 
tonic insulator was first noticed and exploited in a pio- 
neering work of Keldysh [14,15], and explored later by 
others [16,9]. The analogy between different types of de- 
coherence processes acting on polariton condensate, and 
the problem of magnetic and potential impurities in su- 
perconductors suggests that similar methods to that used 
by Abrikosov and Gor'kov (AG) in their theory of gap- 
less superconductivity [1] can be useful in studying the 
properties of polariton condensates in the presence of de- 
coherence. This analogy is however not exact and wc 
will discuss a few essential differences between our the- 
ory and the Abrikosov and Gor'kov approach at the end 
of Section III. An outline of the rest of the paper is as 
follows: In section III we provide details of the method 
in the path integral formulation. In section IV we discuss 
ground state properties of the system in the presence of 
two different types of dephasing while in section V we 
study the excitation spectrum and construct a phase di- 
agram as a function of excitation density, pair-breaking 
decoherence strength and the inhomogeneous broaden- 
ing of energy levels. In section VI we discuss a crossover 
between an isolated condensate and a laser and in sec- 
tion VII the magnitude of the energy gap. In section 
VIII we comment on the applicability of the model and 
the method presented in this work and indicate a few 
directions in which the model could be easily extended. 
Finally, in section IX we discuss recent experiments in 
the light of theoretical predictions presented in this work 
and in section X we briefly summarise the results. 

III. PATH INTEGRAL FORMULATION 

To construct a theory of the coupled system, we will 
exploit an approach based on the coherent state path 
integral. As well as providing access to the mean-field 
equations of the system, such a framework provides the 
potential to explore the influence of fluctuation phenom- 
ena. Working in the grand canonical ensemble, a chem- 
ical potential [i can be used to fix the total number of 
excitations n ex . The quantum partition function of the 
system Z = Tr e -/3( ff -M«ex) can b e expressed as a coher- 
ent state path integral over Fermionic and Bosonic fields, 

Z= f DfaflYlDfcrfj] [] D[c t ,k,c hk ]e~ s . (7) 

3 fc,*=l,2 
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As with the Hamiltonian, it is convenient to separate the 
total action as S = S$ + S$b + Sb where 



f 

S s = SV - / dT Yl fa^o~j 1( < 

J ' 



with S^p = dr-ip(d T + uj c — fi)ip denotes the action of 
the internal electron and photon degrees of freedom of the 
system, while the action for the coupling of the system 
to the bath takes the form 



Ssl 



^C lk (d T +U!ik)c ik 
,2 L k 



Combining the Fermionic fields of the two-level systems 
into a Nambu-like spinor, 

fa = ( h aj ) , 

the bare Green function assumes the matrix form 

Gqj = -d T <r - {ej - p/2)a 3 - g^u- - gip<r + . 

where a denote the Pauli spin matrices (with (T = 1) 
which operate in the (b, a) space (hereafter, loosely re- 
ferred to as the 'particle/hole' space). Finally, we have 
denned the symmetric and antisymmetric 'densities' ac- 
cording to the relation p(r)^ 1 ' 2 ^ = J2j fa( T ) a (o,3)fa( T )- 

Although a theory of the symmetric and antisymmetric 
processes can be developed in concert, for clarity we will 
present a detailed derivation of the action of the 'pair- 
breaking' decoherence processes imposed by r' 1 ). Later, 
in Section III A, we will restore the decoherence processes 
affected by T^l Thus, for now, we will use the following 



abbreviation T 



(i) 



Tfc dropping the 'channel' index. 



Being Gaussian in the Bosonic fields c k , the degrees of 
freedom of the bath can be integrated out leading to an 
effective interaction of the two-level systems which takes 
the form J D[c k , Ck]e~ SsB = e~ s 'sB, where 



S' SB = I drdr'p(T)^lDk(r 
Jo I 



t'Mt'), (8) 



with D^ 1 = —d T — LUk representing the free propagator of 
the environment. Transforming (8) to the Fourier Mat- 
subara frequency representation, and summing over the 
internal degrees of freedom of the bath, 



~ E v l D k{iv n ) = Mi"n), 



(9) 



where D k (u n ) 



-u>k and v n = 2irn/f3, the induced 



S 'sb = fr(iv n )p{iv n )p{-ivn) 



f 

Jo 



dTdT'f T {T - t') ]T Tr^(r) ® ^.(t')^'(t') ® fa(r). 
n' 



In particular, it can be seen explicitly that the interac- 
tion with the environment introduces an effective quartic 
interaction between the different two-level systems. This 
contrasts with the Maxwell-Bloch equations (6) from 
which one can infer only a lifetime for excitations [12]. 

To develop a mean-field theory of the coupled system, 
it is helpful to affect a Hubbard-Stratonovich decoupling 
of the interaction. Introducing the 2x2 component ma- 
trix field Qjj'(r, r'), which inherits the symmetry of the 
dyadic product <j>j{r) ® <pj> (r'), the interaction generated 
by the bath can be decoupled as 

e~ s '^ = J DQe~ Sq 

f drdr'Mr - r') ]T fa {r)Q jr (r, r')^ (r') 

■'0 



x exp 
where 



S Q - f dTdr'Y^Hr - r')Tr Q lf {T,r')Q r » 

with trace taken in the particle-hole space. Combined 
with Ss, an integration over the Fermionic degrees of 
freedom <j> obtains the quantum partition function 



DQe- 



(10) 



interaction assumes the form 



where S = S Q + S,p - Tr In G"\ 

G7/(r,r') = G^5 jr 5{r r') f r (r t')Q 33 ,(t,t'). 

(11) 

and the trace now runs over time and site indicies as well 
as the particle/hole space. At this level, the analysis is 
exact. 

A. Saddle-point equations 

To develop the quantum partition function, further 
progress is possible only within a saddle-point approx- 
imation. Varying the action with respect to Q, one finds 
that each electronic excitation is coupled to the aver- 
age field created by all of the other excitations. In this 
sense, the saddle-point analysis corresponds to a mean- 
field treatment of the interaction between electronic ex- 
citations. This analysis becomes exact when there is a 
large number of electronic excitations coupled to a small 
number of field modes, since the fluctuations of the field 
are then negligible [9]. The mean-field treatment of the 
interaction becomes exact in the thermodynamic limit 
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when N — > oo. The fluctuations above mean-field are of 
the order of 1 / VN. 

In the saddle-point approximation, it is assumed that 
the dominant contribution to the quantum partition 
function (7) arises from those configurations ip and Q 
which minimise the total action. Varying the action S 
with respect to Q one obtains the matrix equation 



Qjj'( T > T ') = 7> G n'{ T , T '), 



(12) 



while varying with respect to ip, the saddle-point equa- 
tion takes the form 



(13) 



Maintaining the analogy with the superconductor, the 
first equation which identifies Q as the self-energy, de- 
scribes the self-consistent Born approximation for the 
Green function, while the second equation represents the 
gap equation for the superconducting order parameter. 
Substituting Eq. (12) into (11) the saddle-point equation 
assumes the form of a Dyson equation 



GTh T , T ') = G7 1 (T,T')S jr -E jr (r,r'), 



(14) 



where 



^(r,r') = -/r(T-T')G jy (T,r') 

denotes the self-energy 

Until now, we have focussed on the impact of the 'pair- 
breaking' perturbation affected by the matrix elements 
Consideration of the symmetric perturbation I^ 2 ' 
follows straightforwardly. In doing so, it may be con- 
firmed that the structure of the saddle-point equations 
are maintained while the self-energy takes the form 

^n> (t, t 1 ) = \f T (r - r> 3 G M v (r, r')a 3 



In steady state one expects the solution of the saddle- 
point equations to depend only on r - r'. In this case, 
a transformation to Matsubara frequencies leads to the 
relation 

= \Y, fr{i<)G{iv n - iv' n ). (15) 



Generally, the solution of the saddle-point equation de- 
pends sensitively on the particular spectrum of decoher- 
ence fr and must be determined self-consistently. How- 
ever, an explicit solution to the saddle-point equation can 
be established in various limits. 

Generally, in order to determine fr from Eq. (9) one 
can assume that coupling constants of the system to the 
bath T k , as well as the bath density of states N(uik), are 
continuous functions of frequency. In this case, one can 



replace the summation over k with an integral over tv k 
(Tk i-> T(u>k), J2k ^ / duo k N k (u k )) whereupon 



fr(w 



iv n ) = j 



dio k 



N k (u k )Y 2 (uj k ) 



-IVr, 



Now, in general fr will exhibit a particular frequency de- 
pendence determined by the density of states N(uj k ) of 
the bath and coupling constant T(uj k ). However two spe- 
cial cases present themselves: Firstly, if we assume that 
both T(uk) and N(w k ) are largely independent of fre- 
quency over a wide range, one finds fr(iv n ) — i2irY 2 N . 
This corresponds to a Markovian approximation to the 
bath in which /r(r) = i2-f 2 S(r). In the second limit, 
if one assumes that the matrix elements and density 
of states are concentrated at zero frequency, one has 
/r(iv„) = 2^ 2 5u n .o- Here, in the static limit, /r(r) = 2^ 2 
is real. Applied to the self-energy, the static limit leads 
to 

Y.jj<(iv n ) = ~f 2 G 33 >(iv n ), 



while, in the Markovian approximation, 
s jy( T > T ) = il 2 G n >(T,r), 



(similarly for the perturbation I^ 2 )). In these two lim- 
iting cases, the saddle-point equations admit a straight- 
forward analytical solution. Here we present a detailed 
analysis for the static limit. 

Noting that photon field in the gap equation (13) cou- 
ples only to the diagonal elements of the Green function, 
Gjj , this suggests a mean-field Ansatz in which the ma- 
trix elements off diagonal in the j space are taken to be 
zero, while the only time dependence of the field is as- 
sociated with oscillation at the chemical potential, fi. In 
this case, the coupled equations (13), (14) can now be 
solved. 

To simplify the algebra we employ a similar mathe- 
matical trick to that used by Abrikosov and Gor'kov in 
their theory of gapless superconductors. Since the over- 
all phase of a coherent state is arbitrary, we can choose 
the mean-field (ip) to be real and present the total Green 
function G (14) in the same form as the zero-order Green 
function Go, 

Gjj(iv n ) = -ivj.n - (ij - fi/2)a 3 - g{i>j)ai, (16) 

using the frequency dependent, renormalised Dj in , ij and 
(ipj). Comparing (14) with (16) we obtain for both type 
1 and type 2 decoherence three equations determining 
the renormalised frequency, energy and coherent photon 
field 



v j,n — v n 7l,2~3 



^ n + fe-M/2) 2 + gW 



e J = £j +7l,2T 



^ + fe-M/2) 2 +.9 2 (^) 2 
9j$j) 



i>l n + (i 3 -^/2) 2 +g 2 (^) 2 



(17) 
(18) 

(19) 
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while the gap equation (13) takes the form 



2(w c - /i) 



(20) 



The average coherent polarisation of the medium can 
be determined from the off-diagonal part of the Green's 
function G (16) 



;• (21) 



Thus, at the mean-field level, we can see from Eqs. (20) 
and (21) that the two order parameters, the coherent 
polarisation and the coherent photon field, are coupled 
according to the relation 



:( p )- 



(22) 
I 



The ratio between the two order parameters is deter- 
mined by the chemical potential, which in the steady 
state can be calculated from Eq. (3). 

The number of electronic excitations, refered to later 
as inversion, can be obtained from the diagonal elements 
of the Green's function 



I(£ (& t 6 ._ a t aj)) 

3 



2(i 3 - n/2) 



"n ,3 3 



^L + fe-M/2) 2 + ff 2 (Vi) 2 



(23) 



Using Eqs. (17)-(19) we can determine the renor- 
malised parameters v^ ni ij and (tpj) as a functions of 
the bare parameters v n , Cj, (tp) and 7 . In the case of 
type 1 decoherence processes we obtain 



< + ^-16 7 2 £ 2 + (£ 2 + 4 7l 2 + ^, 



(24) 



and while for the type 2 decoherence processes we have 



2(^) - W 



(25) 

(26) 
I 



4(^+^)2) 



^E] 2(e - m/2) 2 + 4 7l 2 - vl + ^16 72 2 (^ + S W) + - 4 7l 2 + (27) 



and 



2<V>;> ^ W 



where for both cases 



^ = \/fe-AV2) 2 +ffW 



(28) 
(29) 

(30) 



Substituting Eqs. (24)-(26) or (27)-(29) into Eq. (21), 
summing over the Matsubara frequencies, and using 
Eq. (22) we can determine the ground state coherent 



polarisation (P) and the coherent photon field (ip) as 
functions of the system parameters e, u) c , decoherence 
parameter 7 i j2 and chemical potential \i. The chemical 
potential can be then obtained from Eq. (3). The inte- 
grals over the Matsubara frequencies at zero temperature 
in (21) and (23) as well as the determination of the chem- 
ical potential from (3) have to be performed numerically. 



B. Density of states 

The excitation spectrum, i.e the density of states, can 
be obtained from the diagonal part of the Green func- 
tion on the real frequency axis. Considering the analyt- 
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ical continuation of Gjj{iv n ) — * Gjjiy)i where v n and v 
are the Matsubara and the real frequencies respectively, 
and thus using the usual substitution iv n — ► — v + i5, 
we obtain the relationship between the Green function 
Gjj(iv n ) and the density of states A(v) 

A{v) = y~] lim JiaGjj(-v + iS + fj,), (31) 

which has the following form 

A(v) = Y Im ^ + fe ' M/2) ■ . (32) 

Y ^ + fe-/V2) 2 +.9 2 (^) 2 

Generally, ej and (ipj) are functions of j/, ej, {ip) and 
7 which can be determined from Eqs. (24)-(26) or (27)- 
(29) by the following substitution iv n — > — v + iS. It can 
be shown from Eq. (32) that the system of two-level os- 
cillators with uniform energies, ej = e, in the presence of 
the type 1 processes has a gap, A, in the density of states 
of magnitude 

A = 2V(e-M/2) 2 + flW-47i- (33) 

At very high excitation densities the gap is proportional 
to the coherent field amplitude A w 2g(tp) — 471. At 
very low excitation densities, when {ip) — > 0, we recover 
conventional polaritons for which the chemical potential 
is 2e — (i — > g and the gap A — » g — 471 . 

The major difference between this work and the AG 
theory [1] is that the system studied here has two or- 
der parameters connected through the chemical poten- 
tial which needs to be determined. We use a different 
form of the density of states for the two-level oscilla- 
tors than in their theory. Instead of a flat distribution 
of energies from —00 to +00 used in the AG method 
we first perform the calculations for the degenerate case 
where all two-level oscillators have the same energy e 
and then we use a realistic Gaussian distribution of en- 
ergies, present in the real microcavities. To account for 
these differences we need to include the additional, third 
equation for renormalised e not present in the original 
Abrikosov and Gor'kov method and the constraint equa- 
tion for n ex . In the Abrikosov and Gor'kov theory they 
consider free propagating electrons with momentum k 
over which all the summations are performed. In our 
model of the localised two-level systems the summations 
are performed over the sites where the two-level oscilla- 
tors can be present. Dynamic impurities can easily be 
included in our formalism. 

To perform the calculations we rescale the coherent 
fields by VN and consequently the inversion and the 
number of excitations by N introducing the excitation 
density p ex — n ex /N. In this terminology the minimum 
p ex — — 0.5 corresponds to no photons and no electronic 
excitations in the system. The condition p ex = 0.5 in 
the absence of photons would correspond to all two-level 
oscillators in excited states, thus to the maximum inver- 
sion. 



We calculate first the ground state coherent field 
(■0), the coherent polarisation (P), the inversion and 
the chemical potential as functions of the decoherence 
strength 7 and the excitation density p ex for different 
distributions of excitonic energies. Then we study the 
excitation spectrum of the system for different regimes. 
The ground state properties and the excitation spectrum 
allow us to obtain a phase diagram for different excitation 
densities and decoherence strengths. We consider the in- 
fluence of both type 1 and type 2 decoherence processes 
as well as inhomogeneous broadening of exciton energies. 



IV. THE GROUND STATE — COHERENT 
FIELDS 

A. Type 1 (Pair-Breaking) Decoherence Processes 

To examine the ground state properties of the system 
in the presence of the type 1 decoherence processes we 
study the mean value of the annihilation operator of the 
field and the polarisation. This mean is non-zero only in 
a coherent state. Figure 2 (upper panel) shows the be- 
haviour of the coherent part of the photon field {ip) as the 
decoherence strength 7 is changed for different excitation 
densities p ex and different inhomogeneous broadenings of 
exciton energies. 

For small values of 7/5, up to some critical value 7ci, 
{ip) is practically unchanged while for 7/5 > 7^1 the co- 
herent field is damped quite rapidly with increasing deco- 
herence. This critical value of the decoherence strength, 
7ci is proportional to p exi suggesting that for higher exci- 
tation densities the system is more resistant to dephasing. 
At low excitation densities, where p ex < 0, there is a sec- 
ond critical value of the decoherence strength, 7c 2 , where 
both coherent fields are sharply damped to zero. As the 
excitation density is increased, precisely at p ex — 0, 7^2 
diverges and does not exist for p ex > — coherent fields 
although reduced are never completely suppressed. 

The behaviour of the electronic inversion, which is a 
measure of the excitonic density given by Eq. (23), is 
presented in Fig. 2 (middle panel). In our terminology 
an inversion of -0.5 corresponds to no excitons while an 
inversion equal to zero means that the excitonic system 
is half occupied (roughly 0.5 per Bohr radius in a model 
where double occupation of excitonic sites is not allowed) . 
Inversion of 0.5 would then correspond to 1 exciton per 
Bohr radius. In this region of the decoherence strength 
where {ip) is damped, the inversion increases. At low 
excitation densities (p ex < 0) the inversion approaches 
p ex for 7/g = 7c2 and stays constant as 7 is further 
increased. At high excitation densities (p ex > 0) the in- 
version asymptotically approaches zero with increasing 
dephasing. The maximum value of electronic inversion, 
for any exciton density and decoherence strength, is zero, 
which corresponds to a half filled excitonic system. This 
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is a consequence of our assumption of thermal equilib- 
rium in the exciton-photon system. 



The ratio of coherent polarisation to coherent field 
(P)/(ip) is presented in Fig. 2 (lower panel). For an 
isolated system, where 7 = 0, this ratio depends on the 
excitation density. The condensate becomes more pho- 
ton like as p ex is increased due to the phase space filling 
effect. For finite 7, at a given excitation density, this 
ratio decreases with increasing 7 meaning that the co- 
herent polarisation is more heavily damped than the co- 
herent photon field by the type 1 decoherence processes. 
At p ex < this ratio becomes undefined for j/g > -fc2 
when both coherent fields vanish. 



In order to study the system of realistic, inhomoge- 
neously broadened two-level oscillators we have replaced 
the summations over sites with integrals over the energy 
distribution. We have assumed this distribution to be a 
Gaussian with mean eo and variance a. Our results, pre- 
sented as dashed and dashed-dotted lines in Fig 2 show 
that a Gaussian broadening of energies does not make any 
qualitative difference to the degenerate case. The coher- 
ent fields and the critical values of decoherence strength, 
7ci and -fc2 are, as expected, slightly smaller than in the 
degenerate case but all the regimes are analogous. 




0.5 1 
Decoherence Parameter y/g 



1.5 



FIG. 2. Coherent photon field (if)) (upper panel), inversion 
(middle panel) and ratio between coherent photon field and 
coherent polarisation (P}/(if)) (lower panel) as functions of 
the pair- breaking decoherence strength, y/g for different ex- 
citation densities, p ex and variances of inhomogeneous broad- 
ening a = (solid lines), a — 0.5 (dashed lines) and a — 1.0 
(dotted lines). 
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B. Type 2 (Non-Pair-Breaking) Decoherence 
Processes 




-0.5 1 — ' — ' — ' — ' — 1 — ' — ' — ' — ' — 1 — 1 — ' — ' — L 
2.5 1 — 1 — ■ — ■ — 1 — i — 1 — ■ — ■ — t — i — ■ — 1 — ■ — r 



2- 




)i — i — ■ — ■ — i — i — i — ■ — ■ — i — i — ■ — i — ■ — i — i 
0.5 1 1.5 

Decoherence Parameter y/g 



FIG. 3. Comparison 

between the influence of a pair-breaking (solid line) and a 
non-pair-breaking (dotted line) decoherence processes on (i[>) 
(upper panel), inversion (middle panel) and (P) / '{ip) (lower 
panel) for three different values of p ex . 

We now repeat the analysis for type 2 decoherence pro- 
cesses (5), which act in an exactly the same way on the 



upper and the lower levels of the two-level oscillator. Re- 
call that such processes mirror the effects of non-magnetic 
impurities in the superconductor. 

In Fig. 3 we present for comparison the coherent pho- 
ton field (upper panel), the inversion (middle panel), and 
the ratio between the coherent polarisation and the co- 
herent photon field (lower panel) in the presence of the 
type 1 (solid line) and the type 2 (dotted line) decoher- 
ence processes for three excitation densities p ex = —0.4, 
Pex = — 0-2 and p ex = 0.2. It is evident that the singular 
behaviour at jci an d 7c2 discussed for the type 1 deco- 
herence processes is not present for the case of the type 2 
processes. With an increase of the decoherence strength 
from zero the coherent fields are slightly reduced and they 
slowly decrease, asymptotically approaching zero at low 
excitation densities (p ex < 0) or a constant value at high 
densities (p ex > 0). Although both coherent fields are 
reduced the behaviour of their ratio strongly depends on 
the excitation density. We will show in Section V B that 
the type 2 processes give rise to the broadening of ener- 
gies and then the behaviour of the ratio (P) / (f/>) depends 
on the position of the chemical potential with respect to 
the energy distribution. The type 2 processes can make 
the condensate more photon or more exciton like depend- 
ing on the parameters of the system. The two different 
cases are presented in Fig. 3 (lower panel). 



V. EXCITATION SPECTRUM AND THE PHASE 
DIAGRAM 

To understand the behaviour of the coherent fields we 
study the excitation spectrum of the system in different 
regimes. The density of states of the system with uni- 
form energy distribution for six different values of 7 at 
low excitation density p ex = —0.4 is presented in Fig. 4 
while at high excitation density p ex = —0.2 in Fig. 5. 



A. Type 1 (Pair-Breaking) Decoherence Processes 

In the absence of decoherence (Fig. 4 a and 5 a) we 
have two sharp peaks at two quasi-particle energies, ±E, 
given by equation (30) for the degenerate case (Ej = E) . 
As 7 increases these two peaks broaden, which causes a 
decrease in the magnitude of the energy gap (Fig. 4 b, 
c and 5 b, c, d ). The magnitude of the energy gap in 
the degenerate case is equal to 2E — 47, which is given 
in more detail in equation (33). Finally, precisely at 7ci 
(shown in Fig. 2), these two broadened peaks join to- 
gether and the gap closes (Fig. 4 d and 5 e). When the 
decoherence strength is increased further (Fig. 4 e) these 
two peaks overlap more and the shape of the gapless den- 
sity of states changes. For 7/5 > -fc2 at low densities the 
coherent fields are suppressed, thus Fig 4 f shows the nor- 
mal state density of states in the absence of coherence. 
Finally, at high densities (Fig. 5 f), as one would expect 
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for a laser system, the coherent field is present without a 
sign of a gap. 




y/g = 0.1 
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- y/g = 0.25 
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FIG. 4. Density of states for p ex = —0.4 and different de- 
coherence strengths, 'y/g for a pair-breaking (solid line) and 
a non-pair-breaking (dotted line) decoherence processes. 

Fig 6 shows the density of states for the system of 
the inhomogeneously broadened two-level oscillators with 
standard deviation a = 0.5g for different values of the 
type 1 decoherence strength at two different excitation 
densities. The broadening of the density of states and 
the suppression of the energy gap can be observed as 7 is 
increased. The last curve in the Fig 6 (left panel) shows 
the normal state density where the coherent fields are 
suppressed. In the Fig 6 (right panel) where the coher- 
ent fields are present for all the values of 7 the last curve 
has no sign of a gap. 

There are clearly three different phases depending on 
the decoherence strength, 7 (homogeneous broadening), 
inhomogeneous broadening, a and the excitation density 
p ex . In the Fig. 7 we present a phase diagram for the 
system. The phase boundaries are defined by 701 and 
7C2 for different values of p ex and a. 

Below the white surface, for small decoherence, we 
have a phase in which both coherent fields and an en- 
ergy gap in the density of states are present. In this 
region coherent fields are protected by the energy gap 
and remain practically unchanged while the energy gap 
narrows as the decoherence is increased. At low densi- 
ties within this phase we have a BEC of polaritons with 
the electronic and photonic parts comparable in size. At 



high densities we have a BCS-type of condensate with 
photon component increasing with the excitation den- 
sity. Despite the predominantly photon-like character of 
this phase at very high densities, the coherence of the 
medium is still large and this phase can be distinguished 
from a laser by the presence of a gap in the density of 
states (Fig 6 - right panel). The crossover between high 
and low densities is a smooth evolution and there are no 
rapid phase transitions. 





_ y/g = 1.0 


1 
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FIG. 5. Density of states for p ex = 0.2 and different deco- 
herence strengths, -y/g for a pair-breaking (solid line) and a 
non-pair-breaking (dotted line) decoherence processes. 

Inhomogeneous broadening has some influence on this 
phase mainly at low excitation densities. The critical de- 
coherence 7ci, which defines a boundary between gapped 
and gapless phases, decreases with an increase in the in- 
homogeneous broadening a. At high excitation densities 
inhomogeneous broadening has very weak influence on 
the system. 

Between the white and gray surfaces there is a phase 
where the coherent fields are present without an energy 
gap in the density of states. This phase exists for all 
values of the decoherence parameter at high excitation 
densities p ex > 0. The coherent fields are no longer pro- 
tected by the gap and get reduced as the decoherence 
is increased. The coherent polarisation is more heavily 
damped than the coherent field. Within this phase we 
have at low densities a gapless, light-matter condensate, 
analogous to gapless superconductivity, whilst at very 
high densities we have essentially the laser system. Fi- 
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FIG. 6. Density of states for the Gaussian broaden case with a = 0.5 for different values of the pair-breaking decoherence 
strength at p ex = —0.3 (left panel) and p ex = 0.2 (right panel). 



FIG. 7. Phase diagram. The phase boundaries 701 and 
7C2 as marked in Fig 2. 




B. Type 2 (Non-Pair-Breaking) Decoherence 
Processes 

The density of states in the presence of type 2 pro- 
cesses is shown in Figs 4 and 5 (dotted lines). It can be 
seen that, although the two quasiparticle peaks get very 
broad, the gap is only slightly affected by the type 2 pro- 
cesses and is always present. Even for much larger values 
of 7 than presented in Figs 4 and 5 the gap is not sup- 
pressed. The gap in the density of states is present until 
the coherent fields get completely suppressed. At high 
excitation densities (p ex > 0) coherent fields are always 
present and thus the density of states will have a gap for 
all values of the decoherence strength. 

The type 2 processes give a similar effect as the inho- 
mogeneous broadening of energy levels in the case of an 
isolated system (see Ref. [9]). The difference is that the 
density of states in the presence of the type 2 processes 



has sharp boundaries. This is a result of the method 
which is used to perform the calculations. It can be seen 
that if the bath operators for the type 2 processes, C2's in 
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expression (5), were just numbers, this term could have 
been included into the second term in equation (1) and 
would had just given a random ej. In our calculations the 
C2's are operators and we use a self-consistent Born ap- 
proximation, which is not exact, and does not correctly 
reproduce the tail of the distribution. Thus, because of 
the way the calculations are done, the density of states 
produced always has sharp boundaries. 

In the Abrikosov and Gor'kov theory, due to the 
flat density of states used in the calculations, the non- 
magnetic impurities do not influence the superconducting 
state at all. In the case of degenerate or realistic, Gaus- 
sian distributed energies of the two level oscillators the 
type 2 processes have some quantitative influence on the 
coherent fields and the gap but cannot cause any phase 
transitions. 



VI. CROSSOVER TO LASER 

The features of the second phase (between the white 
and the gray surfaces) in Fig. 7 at high p ex are essentially 
the same as those of the laser system. The laser operates 
in the regime of a very strong decoherence, comparable 
with the light-matter interaction itself. In the presence 
of such a large decoherence laser action can be observed 
only for a sufficiently large excitation density. The coher- 
ent polarisation in a laser system is much more heavily 
damped than the photon field and the gap in the density 
of states is not observed. Thus the laser is a regime of 
our system for very large p ex and 7. 

Laser theories, due to the approximations on which 
they are based, can only be valid in a regime where the 
gap in the density of states is suppressed and thus for a 
large decoherence. At the time when these theories were 
proposed they were deemed sufficient as most lasers oper- 
ate in such a regime. Miniaturisation and improvements 
in the quality of optical cavities in recent years can lead 
to a large suppression of decoherence in laser media. For 
small decoherence and very small pumping in comparison 
to decay processes, when p ex is much smaller then 0, the 
laser theories would predict a lack of coherence while the 
real ground state of the system would be a more matter- 
like condensate. Thus an extension of laser theories to 
account for the gap in the excitation spectrum and co- 
herence in a media is necessary. In our theory the laser 
emerges from the polariton condensate at high densities 
when the gap in the density of states closes for large 
decoherence and thus is analogous to a gapless supercon- 
ductor. 

When both the gap and the coherence in a media are 
taken into account, in contrast to a traditional laser, the 
coherent photon field can be present without the popu- 
lation inversion in the media. The polariton condensate 
is thus an example of a laser without inversion. 

However, it has to be pointed out that there is no for- 
mal distinction between a laser and a Bosc condensate of 



polaritons. In the laser, coherence in the medium (man- 
ifested by the coherent polarisation) although small, is 
not completely suppressed so the laser can be seen as 
a gapless condensate with a more photon-like character. 
One of the possible distinctions between BEC and laser 
could be an existence of an energy gap in the excitation 
spectrum. 



VII. THE GAP 




FIG. 8. Energy gap Ao in the units of g in the absence 
of decoherence as a function of excitation density p ex , and 
Gaussian broadening, a. 

In the degenerate case the magnitude of the gap in 
the density of states is given by (33) If we consider an 
inhomogeneous broadening of excitonic energies with a 
Gaussian distribution then there will always be some os- 
cillators with energies which make the term (e — /u/2) 2 
vanish and thus, strictly, the magnitude of the gap would 
be 2g{ip) — 471 at all excitation densities. This means 
that the gap would vanish at (tp) — > in contrast to de- 
generate case in which its magnitude would be g — 471. 
However, the contribution of these states, arising only 
from the tails of the Gaussian distribution, would be very 
small essentially leading to a soft gap of the magnitude 
of g — 471 at very low densities as in degenerate case. 

We call A the gap in the absence of decoherence when 
71 = 0. Omitting the case when n exc — 0, discussed 
above, we calculate numerically the magnitude of the 
gap, in the units of the dipole coupling g, for a wide 
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range of excitation densities and inhomogeneous broad- 
ening a > 0. A very much depends on the excitation 
density. The Gaussian broadening of energies decreases 
the gap at low densities but has almost no influence at 
high excitation densities. In the presence of the decoher- 
ence (homogeneous broadening) the actual gap A would 
be 

A = A - 4 7l (34) 

Fig. 8 and equation (34) can be used to estimate the mag- 
nitude of the gap in particular experimental systems in 
which the values of broadenings, density and the dipolc 
coupling is known. Some estimates of the magnitude 
of the energy gap in particular materials and conditions 
were reported in a brief report [17]. 

VIII. APPLICABILITY OF THE MODEL AND 
THE METHOD 

As already stated in Section II, the model of two-level 
systems interacting only via cavity photon field gives a re- 
liable description of tightly bound, Frenkel-type excitons 
localised by disorder or bound on impurities, molecular 
excitons in organic materials, atoms in the solid state or 
Josephson junctions arrays in microwave microcavities. 
In the classical limit HmN^oo((nex} /N) = p — > const, it 
has an exact solution of the mean-field form 

\\,u,v) = e^l[(v J bl+u 3 al)\0). (35) 

3 

However, the case of Wannier excitons propagating in 
the sample, present in some clean semiconductors (for 
example GaAs quantum wells), is conceptually not very 
different. The sums over sites in the microscopic Hamil- 
tonian would have to be replaced with sums over mo- 
mentum states and Coulomb interactions between elec- 
trons and holes would have to be included. Such a model 
without photons was studied by Keldysh and Kopaev [14] 
and for coherently driven semiconductor by Schmitt-Rink 
et al. [18]. If we treat the Coulomb interaction on the 
mean-field level then the natural extension of the Keldysh 
mean-field wave-function to account for photons will be 

|*o>- e A ^n[ u fe + ^I«fe]|0) , (36) 
k 

which is of analogous form to that for the two-level sys- 
tems case (35). Then the influence of the various deco- 
herence processes can be included in exactly the same 
way as described in Section II. 

At the level of mean-filed, we do not expect that the 
change of the model from that of two-level oscillators to 
free electron and holes with a Coulomb interaction would 
cause any dramatic differences. It was shown within 
the mean-field techniques of Keldysh [14,15,19] that the 
Coulomb interaction has a pairing effect and leads to a 



formation of a coherent excitonic insulator phase and this 
would only be enhanced by the dipole coupling. 

At high excitation densities (large photon fields) the 
dipole interaction between excitons and photons is the 
dominant interaction. In the case of a driven excitonic 
system Schmitt-Rink at al [18] pointed out that for a 
very large pumping, and thus high excitation densities, 
the Coulomb interaction is just a very minor correction 
to the dominant dipole coupling. They also state that, in 
the absence of the Coulomb interaction, their results for 
propagating electrons and holes are equivalent to those 
obtained for an ensemble of independent two-level oscil- 
lators as optical transitions with different k decouple. 

The Coulomb interaction would be more important at 
low densities. However, in this case, the dominant contri- 
bution arising from the interaction between an electron 
and hole within the same exciton, is taken into account in 
model studied in this work. All other Coulomb interac- 
tions are much weaker and, as with the dipole coupling, 
favour condensation. 

Although we do not expect qualitative differences, it 
would be interesting, especially at low excitation den- 
sities, to perform similar calculations in the momentum 
representation with the Coulomb interaction. This would 
allow one to study the relative influence of the two inde- 
pendent pairing mechanisms, namely the dipole coupling 
and Coulomb interactions, on the condensate. 

Summarising, our method of including and studying 
the decoherence effects is general and can be applied in 
the same way as in this work to propagating excitons as 
well as to coherently driven condensates. The coherent 
photon field in the driven system is an external, fixed, 
parameter and not a self-consistent field satisfying equa- 
tion (20). In contrast the gap in the density of states, 
proportional to the coherent field amplitude, is present 
exactly as for an equilibrium condensate. The model 
can be applied to obtain a qualitative description of the 
physical behaviour even for propagating or weakly bound 
excitons. It would give similar predictions to the model 
based on propagating electrons and holes with Coulomb 
interactions treated within the mean-field approximation. 
It does not, however, include screening and other non 
mean-fields effects. 

If the microscopic origins of the decoherence processes 
are known for a particular experimentally studied sys- 
tem, the theory could be extended to include this in- 
formation. The detailed account of the coupling con- 
stants and the density of states for the environment can 
be easily included within this framework. In this work 
we have assumed a static bath or Markovian (Ohmic) 
bath but in general the bath propagators are frequency 
dependent and thus (tp) would depend on frequency as in 
the Eliashberg theory of gapless superconductors [20,21]. 
The phcnomcnological constants 7 i and 72 can thus be, 
in principle, obtained from this microscopic calculation 
for a particular system. 

The results presented in this work are performed at 
zero temperature for which the summations over the 
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Matsubara frequencies in (20), (21) and (23) become in- 
tegrals and are performed numerically. The extension to 
finite temperatures is straightforward. At finite tempera- 
tures the summations over a discrete v n can be performed 
using the standard techniques [22] and thus finite tem- 
peratures could be easily studied. 

A much more important extension than the finite tem- 
perature case is the problem of non-equilibrium systems. 
For a system with very strong pumping and decay pro- 
cesses the thermal distribution of the relevant quasiparti- 
cles cannot be assumed. The influence of non-equilibrium 
distributions on the polariton condensate will be pub- 
lished elsewhere. 



IX. REMARKS ON EXPERIMENTAL 
REALISATION OF BEC OF POLARITONS 

While stimulated scattering into the lower polariton 
branch (reminiscent of Bosc statistics) was clearly ob- 
served and reported by several groups [23-32], there is 
still no evidence of excitonic coherence, and therefore 
BEC of polaritons. The experiments that have been per- 
formed fall into two categories. In the first, polaritons 
arc pumped coherently, with a laser having an angle of 
incidence chosen at a "magic angle" close to the bottom 
of the lower polariton branch. In the second set of exper- 
iments the pumping is incoherent. The laser pumping is 
tuned to be well above the ground state so that polari- 
tons have to undergo several inelastic scattering events 
before reaching the ground state and thus loosing the 
initial coherence of the pump. 

Although not yet measured, the potential coherence of 
polaritons in the experiments of the first category would 
be inherited from the pump and the behaviour could 
be explained in terms of parametric oscillation. Exper- 
iments of the second category could be potential candi- 
dates for spontaneous BEC of polaritons but evidence for 
it are very sparse. The coherence of the photon field emit- 
ted at polariton frequency has recently been measured in 
one experiment [28] but this alone does not prove that 
there is a coherence in the excitonic part. A laser is an 
example where coherent photon field is generated with- 
out large excitonic coherence. A more direct evidence of 
excitonic coherence is necessary. Large excitonic coher- 
ence would map to a large gap in the excitation spectrum 
which should be seen in the the incoherent luminescence. 

Experiments on cavity polaritons which report the 
stimulated scattering into the lower polariton branch are 
performed at low densities to avoid the fermionic phase 
space filling effect, so that two clear polariton peaks could 
be seen. For such densities the gap in the excitation 
spectrum would be very small and easily suppressed by 
the decoherence processes in the sample. Indeed such a 
gap has never been observed in the photoluminescence 
spectrum. The attempt to increase the density of polari- 
tons in these experiments results not in the formation of 



the condensate but in a switching into the weak-coupling 
regime and lasing. 

This is not surprising since the increase in the density 
of polaritons in these experiments is obtained by increas- 
ing the pumping of excitons. Incoherent pumping is a 
pair-breaking decoherence process and thus the increase 
in the pumping intensity results in the increase of not 
only the density of excitons but also decoherent scat- 
tering. As shown by Eastham and Littlewood [9], the 
fermionic structure of excitons does not prevent conden- 
sation, even at very high densities. Thus, in the current 
experiments, we suggest that it is not a phase-space fill- 
ing effect which lead to a laser as the pump intensity is 
increased, but the increase in the decoherence strength. 
If the polariton condensate is ever to be observed an in- 
crease in density without an increase in decoherence is 
necessary. 

Localised and tightly bound excitons, such as exci- 
tons in disordered quantum wells or Frenkel-type excitons 
in organic compands, seem to be better candidates for 
observing a polariton condensate than the high-quality 
GaAs quantum wells with weakly bound, delocalised ex- 
citons. Static disorder is not a pair-breaking effect and 
would have a weaker influence on the condensate than 
screening and ionisation in the case of delocalised or 
weakly bound excitons. Additionally the polariton split- 
ting reported in organic materials is as large as 80 meV 
[5,6] in comparison to the upper bound of 20 meV in 
GaAs. Another very good candidate for observation of 
effects described in this work would be atoms in solid 
state, glass spheres, dilute atomic gases in which the de- 
phasing is particularly weak in comparison to the dipolc 
coupling. 



X. SUMMARY 

We have studied the equilibrium Bosc condensation of 
cavity polaritons in the presence of decoherence. It has 
been shown [12] that the widely used Langevin equations 
with the constant, frequency independent decay rates are 
not valid for systems with large polarisation in which the 
gap in the density of states is present. In the regime 
of weak decoherence these processes have to be treated 
sclf-consistcntly, in such a way that the gap in the den- 
sity of states is taken into account. We have proposed 
a self-consistent method analogous to the Abrikosov and 
Gor'kov theory of gapless superconductivity [1] which al- 
low us to study all regimes of the system as the decoher- 
ence is changed from zero to large values and for a wide 
range of excitation densities. 

We have found that at small decoherence the polari- 
ton condensate is protected by the energy gap in the 
excitation spectrum. The gap is proportional to the co- 
herent field amplitude and thus the excitation density, 
so the condensate, unlike the excitonic insulator [19], is 
more robust at high densities. This gap gets smaller and 
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eventually is completely suppressed as the decoherence is 
increased. We have shown that there is a regime, anal- 
ogous to the gapless superconductor, when the coherent 
fields are present without an energy gap. This regime, at 
very high excitation densities, has most of the features of 
a photon laser. 

We have studied the influence of two different types 
of processes, both pair-breaking and non-pair-breaking 
ones as well as inhomogeneous broadening of the energy 
levels. We have shown that only the pair-breaking pro- 
cesses can lead to phase transitions. Non-pair-breaking 
processes and inhomogeneous broadening of energies can 
quantitatively change the behaviour of the system but 
cannot prevent condensation. 

We have studied a general form of the interactions, in- 
troducing pair-breaking and non-pair-breaking processes 
and thus our theory is applicable to a wide range of sys- 
tems. For more detailed results, the particular origin of 
the interactions and thus the particular density of states 
for the baths have to be taken into account. 

We have studied the whole phase diagram of the sys- 
tem given by the Hamiltonian (4) for different values 
of the decoherence strength, inhomogeneous broadening 
and the excitation densities and established the crossover 
between an isolated polariton condensate and a photon 
laser as the decoherence strength is increased. 

Our results suggest that, in contrast to traditional 
lasers, coherent light can be generated without a popu- 
lation inversion. This work generalised the existing laser 
theories to include the gap in the excitation spectrum 
caused by the coherence in a media in the low decoher- 
ence regime. 
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